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Abstract 

The leading-particle effect and the factorization property of the scattering am- 
plitude in the impact parameter space are used to study semiclassical production of 
pions in the central region. The mechanism is related to the isospin-uniform solu- 
tion of the nonlinear cr-model coupled to quark degrees of freedom. The multipion 
exchange potential between two quarks is derived. It is shown that the soft chiral 
pion bremsstralung also leads to anomalously large fluctuations in the ratio of neu- 
tral to charged pion.. We show that only direct production of pions in the form of 
an isoscalar coherent pulse without isovector pairs can lead to large neutral-charged 
fluctuations. 
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1 Introduction 



The old puzzle in cosmic-ray observations is the existence of few events charac- 
terized by an anomalously large number of charged pions in comparison with the 
number of neutral pions, the Centauros [1], indicating that there should exist a 
strong long-range correlation between two types of the pions. The negative results 
of the accelerator searches for Centauros at CERN [2,3,4,5] and at Fermilab [6] sug- 
gest that the mechanism for their production is not yet well understood and that 
the production threshold must be larger then 1.8 TeV. 

Such long-range correlations are possible if pions are produced semiclassicaly 
and constrained by global conservation of isospin [7-12]. 

Although the actual dynamical mechanism of the production of a classical 
pion field in the course of a high-energy collision is not known, there exist numerous 
interesting recent theoretical attemps to explain Centauros either as different types 
of isospin fluctuations due to the formation of a disoriented chiral condensate (DCC) 
[13-17], or as multiparticle Bose-Einstein correlations (BEC) [18], or as the formation 
of a strange quark matter (SQM) [19]. Among the most interesting speculations is 
the idea of DCC that localized regions of misaligned chiral vacuum might occur 
during the ultrahigh-energy hadronic and heavy-ion collisions. These regions, if 
produced, would behave as a pion laser, relaxing to the ground state by coherent 
pion emission . It is generally accepted that the fluctuation of the ratio of neutral to 
charged pions of the Centauro type could be a sign of the DCC formation provided 
that a single large domain is formed, containing a large number of low px pions. 
Since the pions formed in the DCC are essentially classical they form a quantum 
superposition of coherent states with different orientation in isospin space. If all the 
pions in the domain are pointing in the same isospin direction and the condensate 
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state is a pure isoscalar then the probability distribution of the ratio of neutral to 
total pions in the domain will be given by [10,13,14,16] 

1 



DCC 



(/) 



2v/J 

where / = n /n is the ratio of the number of 7r 's in the DCC divided by the total 
number of pions. The distribution Poccif) is clearly very distinct from the binomial 
distribution which assumes equal probability for production of n + , 7r_ and 7r pion 
and is of the form 

P B (n ,n) = [_ - - (2) 



which in the limit as n — > oo, n — > oo with / fixed, approaches a delta function at 
/ = 1/3. 

The possibility of observing the DCC type fluctuations critically depends on 
the size and the energy content of the DCC domain. If this domain is of the pion 
size, the effect of DCC will be too small to be observed experimentally. Therefore 
the signals such as the isospin fluctuations, the enhancement of the number of low 
Pt pions, and the suppresion of HBT correlations are all characteristics of a single 
large coherent emission domain [20]. The space-time scenario of the formation and 
decay of the DCC is usually studied within one of the simplified versions of the chiral 
effective Lagrangians, either the linear or nonlinear sigma model [17,21]. However, 
the use of sigma models, be they linear or nonlinear, is only a rough approximation 
to the true situation, because the couplings of pions and sigma to the constituent 
quarks may be large and their effect should not be ignored. 

In the early models [7,8] the coherent production of pions was taken for granted 
and considered as a dominant mechanism. These models also predict strong nega- 
tive correlations between the number of neutral and charged pions. In fact, it was 
observed long ago [9,10] that an uncorrelated jet model involving only pions with the 
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global conservation of isospin gives the same pattern of neutral/charged fluctuations 
as observed in Centauro events. This strong negative neutral-charged correlation is 
a general property of the independent pion emission in which the cluster formation ( 
or the short-range correlation between pions ) is not taken into account [11,12,22,23]. 
What makes these older results interesting is the possibility that such a coherent 
isosinglet pulse of pions could be produced in the high-energy collisions as a result 
of de-excitation of a highly excited region of space, in which the chiral orientation 
of the order parameter - the quark condensate - is different from that of the normal 
vacuum configuration [17,24]. 

In this paper, following the approach of our earlier papers [23], we consider 
in Section 2. the leading-particle effect as a source of a classical pion field in the 
impact parameter space. Since the condition of approximately vanishing isospin of 
the coherent pulse of pions is crucial for the observation of the DCC phenomena, 
pions should be produced in the nearly baryon-free central rapidity region . 

In Section 3. we show how the factorization property of the scattering ampli- 
tude in the impact parameter space may be related to the isospin-uniform solutions 
of the equation of motion of the nonlinear sigma model in the presence of a classical 
source containing constituent quarks [25]. We also study in Section 3. the quan- 
tum case of the nonlinear a model and derive the form of the multipion exchange 
potential between two quarks. 

Our results are summarized in the Section 4. Our general conclusion is that 
within the nonlinear a model the large isospin fluctuations are consequence of singly 
produced pions which are constrained by global conservation of isospin. 

2 Pion production from a classical source 



4 



2.1 The eikonal S matrix with isospin 

At high energies most of the pions are produced in the central region. To isolate the 
central production, we adopt the high-energy longitudinally dominated kinematics, 
with leading particles retaining a large fraction of their incident momenta. We 
assume that the collision energy is large enough so that the central region is free of 
baryons. The energy available for the hadron production is 

Ehad = -V~S — Ei ead i n g (3) 

which at fixed total cm. energy y/s varies from event to event. Using the following 
set of (3n + 2) independent variables s, {<&t, Hi} = Qi, i — 1,2, . . .n, the n-pion 
contribution to the s-channel unitarity can be written as an integral over the relative 
impact parameter b of the two incident leading particles: 

En(s) = / d 2 b ]J dq t | T n (s, M . . . n) | 2 , (4) 

where dq = d 2 qTdy/2(2n) 3 . The normalization is such that 

F n (s) = sa n (s), 

CO 

Vinel(s) = J2 a ™( S )- ( 5 ) 

71=1 

The basic assumption of the independent pion-emission model, in b-space is 
the factorization property of the scattering amplitude T n (s, b; 1 ... n) : 

T n (s,b;l...n) = 2sf(s,b) [J Il^=f^M;gJ, (6) 

where, owing to unitarity, 



<r=+,o ,- v=i v Str- 



and 



n 



f(s,b) f= e -n{s,b) (7) 



[s,b) = J2 [dq\J a (s,b;q)\ 2 (8) 
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denotes the average number of emitted pions at a given impact parameter b. The 
functions | J a (s, b; q) | 2 , after the integration over b, control the shape of the single- 
particle inclusive distribution. A suitable choice of these functions also guarantee 
that the energy and the momentum are conserved on the average during the collision. 
The factorization properties of T n in the form (6) is a consequence of the pion field 
satisfying the following equation of motion 

{n + ml)n(s,b;x)=j(s,b;x), (9) 

where j is a classical source related to J(s, b; q) via the Fourier transform 

J(s, b;q)= J d 4 xe tqx j{s, b; x). (10) 

The standard solution of Eq. (9) is usually given in terms of in- and out-fields 
that are connected by the unitary S- matrix S(s, b; J) : 

?W = S^7li n S = 7T in + jfd, (11) 

where 

n d (x) = J d 4 x'A(x - x'- n)j{s, 6; x'). (12) 

and A = A ret — A a d v . The S matrix following from such a classical source is still 
an operator in the space of pions. The coherent production of pions is described by 
the following S matrix: 

S{s,b;J) = Jd 2 e\e)D{s,b;J)(e\, (13) 

where | e ) represents the isospin-state vector of the leading-particle system. The 
quantity D(s,b] J) is the unitary coherent-state displacement operator defined as 

D(s, 6; J) = exp[ J dqf{s, 6; q)d\q) - h.c], (14) 

where cft(q) denotes the creation operator of the pion field. 
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2.2 Distribution of pions in isospace 



If the isospin of the incoming leading particles is II 3 , then the initial-state vector of 
the pion field is S(s, b; J) \ II 3 ), where | II 3 ) denotes a vacuum state with no pions. 
The n-pion production amplitude is 

iT n (s, b;qi... q n ) = 2s(l'l' 3 ; q 1 ...q n \ S(s, b; J) \ II 3 ), (15) 

where I'I' 3 denotes isostate of the outgoing leading particles. 

If each possible isospin (/', of the outgoing leading particle system is pro- 
duced with equal probability, then we can sum over all (I'I 3 ) using the group theory 
alone to obtain the following probability distribution of producing n + n + , n_n-, and 
n 7r°: 

P Il3 (n + n_n )N Il3 = J2 d 2 bdq x dq 2 ■ ■ ■ dq n \ (I'I' 3 n + n_n \ S(s,b;J) \ II 3 ) | 2 (16) 

TV 

n = n + +n_ + n (17) 

where Nn 3 is the corresponding normalization factor determined by 
E{n}Pii s (n+,n_,n ) = 1. 

This is our basic relation for calculating various pion-multiplicity distributions, 
pion multiplicities, and pion correlations between definite charge combinations. In 
general, the probability -P// 3 (n + n_n ) would depend on (I'I' 3 ) dynamically. The 
final-leading-particles tend to favor the (/', I 3 ) w (/, I 3 ) case. However, if the leading 
particles are colliding nuclei, an almost equal probability for various (/', I 3 ) seems 
reasonable aproximation owing to the large number of possible leading isobars in 
the final state. In the model leading to Eq.(16) the dynamical information which 
restricts the possible values (I'I 3 ) was not taken into account. Nevertheless, even 
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if we assume that (7', I' 3 ) « (7, 7 3 ) the pion's cloud could still have the isospin 
7 W = 0,2,...,27. 

Recent studies of heavy-ion collisions at the partonic level [26] argue that the 
central region is mainly dominated by gluon jets. The valence quarks of the incom- 
ing particles which escape from the interaction region form the outgoing leading 
particle system. Since gluon's isospin is zero, it is very likely that total isospin of 
the produced pions in the central region is also zero. This picture is certainly true 
if the central region is free from valence quarks, the situation expected to appear at 
the extremelly high collision energies. 

If the conservation of isospin is a global property of the colliding system, re- 
stricted only by the relation 

j = /' + 4 (is) 

where I n denotes the isospin of the emitted pion cloud, then J(s, b; q) should be of 
the form 

J(s,b;q) = J(s,b;q)e, (19) 

where e is a fixed unit vector in isospace independent of q. The global conservation 
of isospin thus introduces the long-range correlation between the emitted pions. 

We shall analyze the isospin structure of the model in the so called grey-disk 
approximation in which 

n(s,b) =n(s)0(R(s)-b), (20) 

where n(s) denotes the mean total number of pions produced, R(s) is related to 
the total inelastic cross section, and 6 is the step function. For 7 3 = 7, it is a 
straightforward algebra to calculate the probability of creating n pions of which n 
are neutral pions: 

Pi(n | n) = Pn(n + n_n ) (21) 

n++n_=n— no 
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( n \ B(n + ln-n + I + l) 
~ \no) B&I + l) ■ {22) 

r(x)r(v) 

Here B(x,y) = y (x + y) * S ^ e ^ u ^ er beta function. Note that Pi (no \ n) differs 

n 

considerably from the binomial distribution given here by 3 _n 2 n_n °(no). 

In Fig. 1 we show the behavior of Pi for n = 50 and isospin I = 0, 1 of 
the initial-leading-particle system in comparison with the corresponding binomial 
distribution. It is easy to see that in the limit n — > oo, with / = fixed, the 
probability distribution nPi(n \ n) scales to the limiting behavior: 

This limiting probability distribution is different from the usual Gaussian random 
distribution for which one expects Pi(f) to be peaked at / = | as n — > oo. How- 
ever, if pions are produced through the coherent production of clusters which sub- 
sequently decay into two or more pions the fraction of neutral pions in an event 
changes substantially. 

2.3 Production of isovector clusters 

Let us assume that pions are produced both singly and through isovector clusters of 
the p type [23]. In this case, the most appropriate tool for studying pion correlations 
is the generating function Gu 3 (z,n_): 

G Il3 (z,n_)= ]T P Il3 (n + ,n_,n )z no , (24) 

no,n + 

from which we can calculate, for example 

(n ) n _ = ^lnG J73 (l,n_), (25) 

f° n _ = ^\nG Il3 (l,n_), (26) 
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and 

^ = is^ G,lM (27 > 

The form of this generating function in the case of / = I3 = is particularly simple: 

^(^jJ'^-BM), (28) 
Jo n l 

where 

2A(z, x) = (1 - x 2 )^ + ^(l - x 2 )n p + 2x 2 n p (29) 

and 

2B(z, x) = 7^(1 + x 2 - 22x 2 ) + n p (2 - z(l - x 2 )). (30) 

Here n n denotes the average number of singly produced pions, and n p denotes the 
average number of p-type clusters which decay into two short-range correlated pions. 
Note that A(l,x) = B(l,x). 

The total number of emitted pions is 

n = n n + 2n p . (31) 

In Fig. 2. we show the behavior of P(n ) = P o(^o) f° r n = 50 and different 
combinations of (n n ,n p ). We see that Centauro-type behavior is obtained only for 
n n 7^ and n p = 0. Recent estimate of the ratio of p- mesons to pions, at accelerator 
energies, is n p = 0.107?^. The behaviour of (n ) n _ and can be found in [23]. 

2.4 Random source 

In this section we want to study the multiplicity distribution of pions emitted from 
a classical random source. In previous sections we have averaged our results only 
over a chiral orientations of the classical source but not the overall shape of the 
source. We assume that incident leading particles, which can neither be created nor 
destroyed, are acting as a classical random source for pions. 
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The unitary S matrix following from such a classical random source is S(s,b; J). The 
initial-state vector for the pion field is S(s,b; J) \ II 3 ).ln practice, we rarely have 
any information about this initial state. It means that physical quantities should be 
averaged over the choice of source functions present in the initial state. In quantum 
statistics, the ensemble average is usually performed using the density operator, 
which, in our case, is of the form [27] 

p Ih (s,b) = [S(s,b;J) | II 3 )(II 3 | &(s, b; J)] av (32) 

and normalized to unity: Tr[pu s ] = 1 In terms of the pion-number operator 

N = N + + AL + iVo 

= S / d Q a l(Q) a a(q) (33) 
the pion-generating function in the impact-parameter space becomes 

G Ih (z ,z + ,z.) = Tr[p Ih (s,b) : e^ z ^ N ° :] (34) 
= (exp[^2(z a - l)n a (s,b)]) Il3 (35) 

a 

where 

i a (s,b) = J dq | J a (s, b;q) | 2 . (36) 



The global conservation of isospin and the identity of pions require 

J a (s,b\q) = J(s,b;q)e (37) 

where 

e = e 3 (38) 
e± = -j={e 1 ^ie 2 ). (39) 



Since 



n a (s, b) — \ e a \ n(s, b) (40) 
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we can find the generating function in 6-space for producing n c h charged and n 
neutral pions which is of the form 

G Ih (z , z ch ) = (I+ ^j^fyj £ dx \ Pf*(x) | 2 ( e «W>(*>) (41) 

where 

z(x) = x 2 z + (1 - x 2 )z ch - 1 (42) 

and P/ 3 is the associate Legendre polynomial. The case I = I3 = is particularly 
simple. The corresponding generating function is 




where G = Goo and 

*(/) = fzo + (1 - f)zch - 1. (44) 
For simplicity of writing we have left out the (s, 6) dependence in n(s, b). 

If the classical random source is such that its fluctuations about the mean are 
Gaussian[27], the generating function becomes 

G(,„,, cfe ) = {^e^^^ A2 « (45) 

where 

d 2 = (n 2 (s,b))-(n(s,b)} 2 - (46) 

The form of this generating function is such that the Centauro type behaviour of 
P(n ) can be expected only in the limit of a very small dispersion d in the number 
of produced pions in b-space. This problem certainly requires further study, in 
particular, if we go beyond the Gaussian approximation for the fluctuations of the 
source about the mean. 
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3 Relation to the nonlinear sigma model 



In this section we study the relationship that can be established between the non- 
linear sigma model coupled to quarks and our eikonal model with factorization. 

3.1 Quark sources of the pion field 

As is well known, the Lagrangian for QCD with two light up and down quarks has 
an approximate global SU(2)l x SU(2)r symmetry, which at low temperatures, 
is spontaneously broken to SU(2)y by a nonzero value of the quark condensate 
{QlQr), which is regarded as an order parameter of the system. This order parameter 
can be represented as a four-component vector = (a, n) buildt from the quark 
densities. The chiral symmetry then corresponds to 0(4) rotations in internal space. 

The true vacuum of the theory is defined as ((f)) = ((<t),0), with (a) ^ 0. 
In QCD the spontaneous symmetry breakdown leads to nearly massless Goldstone 
bosons (the pions) and gives the constituent-quark mass. At low energies and large 
distances ( momentum scale smaller than lGeV ) the dynamics of QCD is described 
by an effective Lagrangian containing the a, n fields and constituent quarks. Re- 
cently, the idea that extended regions of the misaligned chiral field <fi may be formed 
in a very high energy hadron-hadron, hadron-nucleus or nucleus-nucleus collisions, 
has attracted a lot of interest both theoretically and experimentally [17]. If a single 
large DCC domain forms, its decay into pions can lead to anomalous fluctuations 
in the ratio of neutral to charged pions similar to ones observed in Centauro-type 
events. 

In the DCC dynamics we distinguish three stages: formation, evolution and 
decay stage. In the conventional approach [16] one starts with a chirally symmetric 
phase at T > T c and DCC formation happens as T, due to a rapid expansion or 
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cooling, drops below T c spontaneously breaking the chiral symmetry. 

The evolutionary stage of the DCC is usually described by the classical chiral 
dynamics based on the a-model, mostly the linear a- model. For the purpose of 
comparison with our eikonal model we consider the nonlinear a-model coupled to 
quarks at zero temperature [25] which is expected to describe the late stage of the 
DCC evolution. 

The Lagrangian for the nonlinear a-model coupled to quarks is 

L = &Tr(d^d»U) + q(nd)q - gf^qUq (47) 



where 



U = exp(t l5 ^) (48) 



We parametrize the pion field in the following form 

tt(x) = Un{x)6{x) (49) 

where n(x) is an unit vector which determines the isospin orientation of the pion 
field, obeying n 2 — 1. 

The Euler-Lagrange equations of motion for 9 and n are 
oe -sine cos Qd^n-d^n = -i^n ■ {Qr lb Q) 

d^(sm 2 6n x d^n) = -iy^n x (Qr l5 Q) sin 6 (50) 

(51) 

where Q denotes the constituent quark defined by 

Q = e^q (52) 
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and itlq = gf„ is the constituent quark mass. We treat 

- i^QfiftQ = fa) (53) 

Jit 

as a given classical external source and consider the class of solutions that can be 
rotated into a uniform one, n(x) = constant, known as the Anselm-class of solutions 
[13]. The solutions with constant n can be realized if the source points to a certain 
fixed direction e in the isospace: 

j(x)=j(x)e (54) 
The source term in the equation of motion for n disappears if we choose 

n(x) = e. (55) 
Then the equation of motion for the pion field reduces to 

ne(x)=j(x) (56) 

with 

n(x) = fj(x)e (57) 

We see that this is exactly the equation of motion for the pion field that we 
have used in our eikonal model in order to predict the appearance of the Centauro 
type behaviour in heavy-ion and hadron collisions. This relationship therefore offers 
the possibility for studying the importance of various quark sources responsible for 
the DCC formation. 

3.2 Quantum nonlinear sigma model 

In quantum chiral field theory of the nonlinear a-model the role of a strong coupling 
of (^-interaction has the quantity uiq/ f n . However, in phenomenological analyses 
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the strong coupling constant of qn- interaction is qq. It is connected with toq/ /, 
through the Goldberger-Treiman relation 



where qa denotes the axial vector current constant whose value, if different from 
unity, should come, according to Lehmann [28], from higher orders of perturbation 
theory. Since the Lagrangian of the nonlinear a-model is nonpolynomial, a suitable 
renormalization procedure and operator normal ordering should be formulated [29] . 
In this respect, as a first step, we derive the form of the multipion propagator 
between two quarks. Let 



denote the chiral super field of the pion. The Lagrangian describing the interaction 
of this super chiral field with quarks is 




(58) 



$ = : e i76 ^ - 1 : 



(59) 



L = — 




(60) 



The chiral super propagator of the field <3> is defined by 



A* (x) = <T($Or)$(0))> 



(61) 



Its form, after number of algebraic manipulations, is 



A* (x) = {1 <g> 1 + ^( 75 r) <g) (^r)d A }d A [Ach(— )] + 1 



(62) 



where 



<T(7r i ( a ;)7r i (0))> = ^A( a ;). 



(63) 



and 




1 1 



4-7T 2 x 2 — ie 
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The multipion exchange potential between two quarks is related to the Fourier trans- 
form of A$(rr) in the following way: 

u(p' l )u(p' 2 )K^{q)u(p 1 )u(p 2 ) = u' l u' 2 %{q)u 1 u 2 (64) 

where 



u(p) 



tp + m [ u 



6p 



and 



, ap 

CO = UJ 



e v + m 



with uj*ijj = 1. 

We should also mentioned that the soft chiral pion bremsstralung [30-33], in 
which every incoming and outgoing quark line is replaced by 

— * — * 

q — > exp(i~f 5 — )q (65) 

also leads to the distribution of neutral pion of the form which is typical 

for coherent pion production without invoking the notion of DCC formation. We 
conclude therefore that nP(n \ n) ~ \Jn/no is not always a definite signature of 
DCC formation. 

4 Conclusion 

The result of our analysis are the following: 
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• Within the framework of an unitary eikonal model with factorization and 
global conservation of isospin the Centauro-type behaviour can only be ex- 
pected for pions which are produced singly ( see Fig.l and 2 ). 

• The Centauro-type effect depends on isospin of the initial-leading- particle 
system (see Fig.l). 

• The coherent production of p-type clusters supresses the Centauro-type be- 
haviour (see Fig. 2). 

• A classical random pion source, which has Gaussian fluctuations about the 
mean, also supresses the Centauro-type behaviour. 

• The leading-particle effect and the factorization property of the scattering 
amplitude in the impact parameter space may be related to the isospin-uniform 
solutions of the nonlinear cr-model coupled to quarks. 

• The multipion exchange potential between two quarks is derived which is dif- 
ferent from the one-gluon exchange potential. 

• The soft chiral pion bremsstralung also leads to anomalously large Centauro- 
type ratio of neutral to charged pions. 

Acknowledgment This work was supported by the Ministry of Science 
of Croatia under Contract No. 0098104. 
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Figure captions : 

Fig. 1. The distribution of neutral pions as a function of the isospin of the 
incoming leading-particle system in comparison with the binomial distribution. 

Fig. 2. The curves represent P(no) for different combinations of (n n ,n p ), the 
average number of singly produced pions and the average number of p-type clusters, 
respectively. 
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